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Abstract 

It is found that noncommutative ^7(1) gauge field on the fuzzy sphere is equivalent 
in the quantum theory to a commutative 2— dimensional U {N) gauge field on a lattice with 
two plaquettes in the axial gauge Ai = 0. This quantum equivalence holds in the fuzzy 
sphere-weak coupling phase in the limit of infinite mass of the scalar normal component 
of the gauge field. The doubling of plaquettes is a natural consequence of the model 
and it is reminiscent of the usual doubling of points in Connes standard model. In the 
continuum large N limit the plaquette variable W approaches the identity I2N and as a 
consequence the model reduces to a simple matrix model which can be easily solved. We 
compute the one-plaquette critical point and show that it agrees with the observed value 
q;=k = 3.35. We compute the quantum effective potential and the specific heat for U{1) 
gauge field on the fuzzy sphere Sj^ in the expansion using this one-plaquette model. 
In particular the specific heat per one degree of freedom was found to be equal to 1 in 
the fuzzy sphere-weak coupling phase of the gauge field which agrees with the observed 
value 1 seen in Monte Carlo simulation. This value of 1 comes precisely because we have 
2 plaquettes approximating the NC U{1) gauge field on the fuzzy sphere. 
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1 Introduction 

Quantum noncommutative ( NC ) gauge theory is essentially unknown beyond one-loop 
[1]. In the one-loop approximation of the quantum theory we know for example that gauge 
models on the Moyal-Weyl spaces are renormalizable [2]. These models were also shown to 
behave in a variety of novel ways as compared with their commutative counterparts. There 
are potential problems with unitarity and causality when time is noncommuting, and most 
notably we mention the notorious UV-IR mixing phenomena which is a generic property of all 
quantum field theories on Moyal-Weyl spaces and on noncommutative spaces in general [1,3]. 
However a non-perturbative study of pure two dimensional noncommutative gauge theory was 
then performed in [5]. For scalar field theory on the Moyal-Weyl space some interesting non- 
perturbative results using theoretical and Monte Carlo methods were obtained for example 
in [6]. An extensive list of references on these issues can be found in [1] and also in [4] 

The fuzzy sphere ( and any fuzzy space in general ) is designed for the study of gauge 
theories in the non-perturbative regime using Monte-Carlo simulations. This is the point of 
view advocated in [7]. See also [8-10] for quantum gravity, string theory or other different 
motivations. These fuzzy spaces consist in replacing continous manifolds by matrix algebras 
and as a consequence the resulting field theory will only have a finite number of degrees of 
freedom. The claim is that this method has the advantage -in contrast with lattice- of preserving 
all continous symmetries of the original action at least at the classical level. This proposal was 
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applied to the scalar 0^ model in [11] and to the U{1) gauge field in [12] with very interesting 
non-perturbative results. Quantum field theory on fuzzy spaces was also studied perturbatively 
quite extensively. Sec for example [13-15]. For some other non-perturbative ( theoretical or 
Monte Carlo ) treatement of these field theories see [16,24]. 

Another motivation for considering the fuzzy sphere is the following. The Moyal-Weyl 
NC space is an infinite dimensional matrix model and not a continuum manifold and as a 
consequence it should be regularized by a finite dimensional matrix model. In 2 dimensions 
the most natural candidate is the fuzzy sphere Sj^ which is a finite dimensional matrix model 
which reduces to the NC plane in some appropriate large N fiattening limit. This limit was 
investigated quantum mechanically in [14,17]. In 4— dimensions we should instead consider 
Cartesian products of the fuzzy sphere S% [15], fuzzy CP^ [18] or fuzzy S"^ [19]. It is fair 
to mention here that an alternative way of regularizing gauge theories on the Moyal-Weyl NC 
space is based on the matrix model formulation of the twisted Eguchi-Kawai model. See for 
example [20,21,32]. 

The goal of this article and others [12,22] is to find the phase structure ( i.e map the different 
regions of the phase diagram ) of noncommutative U{1) gauge theories in 2 dimensions on the 
fuzzy sphere S%. We consider the fuzzy sphere since it is the most suited two dimensional space 
for numerical simulation because of the obvious fact that it is a well defined object. 

There are three phases of U{1) gauge theory on S^. In the matrix phase the fuzzy sphere 
vacuum collapses under quantum fluctuations and we have no underlying sphere in the con- 
tinuum large N limit. This phenomena was first observed in Monte Carlo simulation in [23] 
and then in [12]. In [22] it was shown that the fuzzy sphere vacuum becomes more stable 
as the mass of the scalar normal component of the gauge field increases. Hence this vacuum 
becomes completely stable when this normal scalar field is projected out from the model. This 
is confirmed in Monte Carlo simulation in [12]. 

In the other phase, the so-called fuzzy sphere phase, there are in fact two distinct regions 
in the phase diagram corresponding to the weak and strong coupling phases of the gauge field. 
The boundary between these two regions is demarcated by the usual third order one-plaquette 
phase transition [25]. This is precisely what we observe in our Monte Carlo simulation of the 
model with a very large mass of the normal scalar field [12]. This result indicates that quantum 
noncommutative gauge theory is essentially equivalent ( at least in this fuzzy sphere phase ) to 
( some ) commutative gauge theory not necessarily of the same rank. This prediction goes also 
in line with the powerful classical concept of Morita equivalence between noncommutative and 
commutative gauge theories on the torus [1,21]. 

In this paper we will give a theoretical proof that quantum noncommutative gauge theory is 
equivalent to quantum commutative gauge theory in the fuzzy sphere-weak coupling phase in 
the limit of infinite mass of the normal scalar component of the gauge field. More precisely we 
will show that the partition function of a f/(l) gauge field on the fuzzy sphere is proportional 
to the partition function of a generalized 2— dimensional U{N) gauge theory in the axial gauge 

= on a lattice with two plaquettes. This doubling of plaquettes is reminiscent of the usual 
doubling of points in Connes standard model [27]. This construction is based on the original 
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fuzzy one-plaquette model due to [26] . 

However in the present article we will show that in order to maintain gauge invariance 
and obtain sensible answers we will need to introduce two different U{1) gauge fields on the 
fuzzy sphere which will only coincide in the continuum large limit. This doubling of fields 
is not related to the above doubling of plaquettes since it disappears in the continuum limit 
where the path integral is dominated by the configuration in which the two U{1) gauge fields are 
equal. Furthermore we will need in the present article to write down two different one-plaquette 
actions on the fuzzy sphere. Linear and quadratic terms in the plaquette variable W are in fact 
needed in order to have convergence of the path integral. We will show explicitly the classical 
continuum limit of these one-plaquette actions. 

Quantum mechanically since the plaquette variable W is small in the sense we will explain 
we can show that the model in the large N limit will reduce to a simple matrix model and 
as a consequence can be easily solved. We compute the critical point and show that it agrees 
with the observed value. We will also compute the quantum effective potential for U{1) gauge 
field on the fuzzy sphere Sj^ in the 1/N expansion using this one-plaquette model. This is 
in contrast with the calculation of the effective potential in the limit — > oo in the one- 
loop approximation done in [22]. The difference between the two cases lies in the quantum 
logarithmic potential which is in absolute value larger by a factor of 4 in the 1/N expansion as 
compared to the one-loop theory. We will discuss the implication of this to the critical point 
and possible interpretation of this result. We will also compute the specific heat and find it 
equal to 1 in the fuzzy sphere-weak coupling phase of the gauge field which agrees with the 
observed value 1 seen in Monte Carlo simulation. The value 1 comes precisely because we have 
two plaquettes which approximate the noncommutative U{1) gauge field on the fuzzy sphere. 

This paper is organized as follows. In section 2 we will briefly comment on the classical 
Morita equivalence between noncommutative gauge theories and commutative gauge theories 
on the torus. In section 3 we will rederive the one-loop result of [22] using an RG method. 
Thus we will explicitly establish gauge invariance of the S|r-to-matrix critical point. Section 4 
contains the main original results of this article discussed in the previous three paragraphs. In 
section 5 we conclude with a summary and some general remarks. 

2 The NC torus and Morita equivalence 

The strongest argument concerning the equivalence between classical noncommutative gauge 
theories and classical commutative gauge theories comes from considerations involving the 
noncommutative torus and Morita equivalence. In this section we will briefly review this result 
following the notations of [1,21]. 

Any U (N) gauge model on the noncommutative torus Tq with a non-zero magnetic flux Q 
can be shown to be Morita equivalent to a C/ (A'o) gauge model on the noncommutative torus 
T3 with zero magnetic flux. The noncommutativity parameter ©' = 2t^9' detT! is given in terms 
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of = 27r9det'E by the equation 



e' = (r 



The integer I which is the ratio / = N/Nq is the dimension of the irreducible representation of 
the Weyl-'t Hooft algebra 



found in the non-trivial solution 
of the model 

Sym 



TiFa = e^"^''T2Ti (2) 



=6^""^^'®!, (3) 



^ J d'xtr^iF,^ - (4) 



Ta are constant SU{N) matrices while aai is a 2 x 2 real matrix which represents the U{1) 
factor of the U{N) group. 

The N X N star-unitary transition functions fla are global large gauge transformations 
whereas fij is a constant curvature on T| which is equal to the curvature of the U{1) background 
gauge field given by 

ai = -^DijX^ , D = 2a^^-^. (5) 

In above we have q/N = m/l where q = Q12 is the component of the antisymmetric matrix 
Q of the non-abelian SU (N) 't Hooft flux across the different non-contractible 2— cycles of the 
noncommutative torus. By construction q is quantized, i.e q E Z. Furthermore m is defined 
by m = |, a; = gcd(g, A^) where 'gcd" stands for the great common divisor. Since I and m are 
relatively prime there exists two integers a and h such that al + hm = 1. These are the same 
integers a and b which appear in (P). 

The period matrix E' of the dual torus T^, is related to the period matrix S of the torus Tq 

by 

S' = -(em + /)S. (6) 

The dual metric is therefore rj' = S's'-^ = (Gm + /)^1 which is to be compared with the original 
metric 77 = SS^ = 1. 

The dual action is by the very definition of Morita equivalence equal to the U{N) gauge 
action Sym on T|, viz 

Sym = -7^ / d^x'trNo{J^iAx'))l (7) 
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where g"^ is given in terms of by the equation 



^'^ = /^(em + /f. (8) 

In other words Sym can also be interpreted as a [/ (A^o) gauge action on T'^, . It is understood that 
the star product here is the one associated with the parameter 6' . It is the U{1) background 
gauge field which is used to twist the boundary conditions on the U{N) gauge field and 
hence obtain a non-trivial field configuration. Indeed the curvature /12 of the vacuum gauge 
configuration Oj on Tg is related to the SU (N) 't Hooft magnetic fiux q by the equation 

If we turn this equivalence upside down then we can obtain a correspondence between a U{1) 
gauge model on a (finite dimensional) fuzzy torus T] and an ordinary U{N) gauge model on 
T^. In particular we remark that if we set 6 = 0, E = 1 and iVo = 1 in the above equations 
then I — N and m — q and we will have an ordinary U{N) on a square torus with non- 
zero magnetic fiux q — and a coupling constant g'^. The dual torus in this case is also 
square since its period matrix is given by — —N.l whereas its noncommutativity parameter 
becomes 

The commutation relation of the NC torus TB becomes therefore 

Z2z[ = z[z2 exp{-27ri^). (11) 

Since the noncommutativity parameter here is rational we know that this Lie algebra must 
have a finite-dimensional NxN representation which can be written down in terms of shift and 
clock matrices as usual. In other words 

4 = , 4 = {Wn)'. (12) 

Vn and Wn are the canonical SU (N) clock and shift matrices which satisfy VnWn = c^WnVn- 
This is indeed a fuzzy torus, i.e T^, = T^. The coupling constant of the U{1) model on is 

9' = g'N. 



3 The Fuzzy Sphere S% 

Let , a = 1, 2, 3, be three NxN hermitian matrices and let us consider the action 

S[Xa] = -jTr[Xa, Xb]' + '-^eabcTr[Xa, Xb]X, + ^TrX^ + MTr{Xlf (13) 
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This action is invariant under U{N) unitary transformations Xa — >UXJJ^ . The trace is 
normahzed such that Trl = A^. a, (3 and M are the parameters of the model. This action 
is bounded from below for all strictly positive values of M. For M = this model is also 
symmetric under global translations Xa — >Xa + Xa'^N where x is any constant vector. We can 
fix this symmetry by choosing the matrices Xa to be traceless. 
The classical equations of motion read 

Ja\X] = Na'lX,, iXab] + 2l3Xa + 2M{Xl X4=0 , a'Xab = X,] + aeabcX,. (14) 

Absolute minima of the action are explicitly given by the fuzzy sphere solutions 

Xa = RLa®ln- (15) 

La are the generators of SU{2) in the irreducible representation -j. They satisfy [La^Li] = 
ieabcLc, Ll = ^{^ + 1)=C2 and they are of size (L + l)x(L + 1), viz N = n{L + 1). R is the 
radius of the sphere given explicitly by the solution of the equation 

(l + ?fl)fl^-afl + | = 0. (16) 

In particular for /5 = M = we have the solution R = a. If we insist that R = a then we will 
have the constraint 

P = -2Ma^C2 = -d?m^ , a = a^N , = (17) 

In general we can show that a solution of (fTHjl exists if and only if (3 is such that /3<4^y+^;?)- 
Explicitly we have R = R\fN = R{a,rn?) with 

5.~ 2^ « + y/g^ - 4(1 + m^)/3 

R^^^^) = ^^(YT^) • (18) 

In the following we will strictly work with the case R= a. We expand around the solution (JT 
by writing 

Xa = aiLa + Aa). (19) 

Aa , a = 1,2,3 are NxN hermitian matrices which admit the interpretation of being the 
components of a U{n) gauge field on a fuzzy sphere of size (L + l)x(L + 1). To see this we 
introduce the curvature tensor by 

a^Fab = i[Xa, Xb] + aeabcXc^Fab = (iCaAb - iCbAa + eabcAc + l[Aa, Ab] ) . (20) 



We also introduce the normal component of A by 

J VV'Slr I , u,— LL I J 

^ 2 V Joi 



= ^L^^$ = \( ^^Aa + AaXa + ^ ] (21) 



where La = [La,--] and the derivations and coordinate-operators on the fuzzy 



sphere Sj^. We can then check that the action ^[Xa] takes the form 

2a W„ ,o 1 



S[Aa] 



AN' 



TrF^b - j^eabcTr[FatA, - -[Aa,A,]A, 



+ -——Tr^^ - -a% - -a%m^- (22) 
iV o 2 



We note that a natural definition of the U (n) gauge couphng constant is given by g"^ 



Also 



we note that Sq = S[Aa = 0] = —\a^C2 — \a^C2m?'- Finally we remark that there is no linear 



term in $. For completeness we will include a linear term in $ as follows 



S[Aa] 



JTrTrFl - eahcTT[FahA, - -[Aa,A,]A 



AN' 
a" 



AN' 

+ ^0. 



:TrF, 



AN 



+ 



ab 



Tr[FabA,--[Aa,At]A 



N 



So 



Aa'^m? 
N 



N 



(23) 



So — Sq 

by 00- 



2a m 0Q. Now in the limit m — >oo the field $ will be equal to a constant given 



3.1 The effective action from an RG method 

We are interested in the partition function 

Z = J [dXje-^[^"] (24) 

For simplicity we consider U{1) theory so that N = L + 1 and the full fuzzy U{1) symmetry is 
given by the gauge group U (N) . The treatement of U (n) is identical. We will also confine our 
analysis to the case where the coupling constant (3 is related to m by (fTTj) . 

We will fix the U{N) symmetry by diagonalizing the third matrix X3. This will clearly 
reduce the original U{N) symmetry group to its maximal abelian subgroup f/(l)^. Although 
this method is not manifestly SU (2)— covariant it is completely gauge invariant since it does not 
require any extra parameter to be introduced in the model unlike other gauge-fixing procedures. 
Thus we will choose a unitary matrix U such that U^X^U = A3 is a diagonal matrix with 
eigenvalues Xa, A = 1,N. We will have the simultaneous rotations U^XiU = Aj, i = 1,2. As 
it turns out X3 = UA^U'^ can also be thought of as a parametrization of the matrix X3 in 
terms of its radial degrees of freedom encoded in A3 and its angular degrees of freedom given 
hj U = e*®. Indeed we can compute the following metric and measure 

Tr{dX,Y = J2d\l + 2j2i>^A- \BfdQABdQ\B 

A A<B 

[dx,] = ( n ^^^) f n - {x{^^/^abdQl^ . (25) 

Va=1 ^ ^A<B ^ ^ ^ 
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The partition function becomes ( since the integration over the unitary matrix U decouples ) 

Z= [[dA,](lld\X-'^^"'-l (26) 



where the action is now given by 

A<B 

SN[Aa\ = -jTrN[Aa,Ab? + '-^eabcTr4Aa,At\A, + PTr^Al + MTrN{Alf. (27) 

We are using now the new notation Tr=TrN. We stress again the fact that this action is 
still symmetric under the abelian U{1)^ transformation Aj — >V~^AiV and A3 — ^A^ where V 
is given explicitly by 

Vab = e''HAB. (28) 

Now we adopt the RG prescription of [28] to find the quantum corrections of this action at 
one-loop. To this end we parametrize the NxN matrices A^ in terms of (A^ — l)x(A^ — 1) 
matrices Da, {N — 1)— dimensional vectors Va and 1— dimensional vectors pa as follows 

Since A3 is diagonal we must have f 3 = while pa = X^. This method consists in finding 
quantum corrections to the action coming from integrating out the 4(A^ — 1) + 3 degrees of 
freedom Vi and pa which we can naturally think of as fluctuations around a background defined 
by the matrices Da- Furthermore it is not difficult to argue that this method is also equivalent 
to the usual Wilson procedure of integrating out the top modes with spin L = — 1 from the 
theory. 

To see this more exphcitly we write {Aa)AB = {Da)AB, {Ai)AN = vf, {Ai)NA = (vf)* and 
{Aa)NN = Pa whcrc A, B = 1, ...,N — 1. We check that the abelian transformations will 
act on Da, Vi and pa as follows 

Da-^WDaW+ , Va-^WVa , -^vtW+ , Pa-^Pa (30) 

where 

{Wa)AB = e-'^'--'-^6AB = e''HV^)AB. (31) 
Next we will denote the (A^ — 1)— dimensional trace by TrTv-i and compute 

Tr^[A,,A,]2 = Trr,.i[Da,D,]^ + vt[-mDj+8D,D,-AD%]v,+Oi3). (32) 

ieabcTrN[Aa, Ab]Ac = ieabcTrN-i[Da, Db]Dc - Qieij:>,vf D^^Vj + 0(3). (33) 
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Tr^Kl = Tr^-,Dl + 2vtv, + pi (34) 

and 

3 

Tr^iAlf = Tr^.,{Dlf + J^p^ + vf [2D%, + 2D,D,]v, + 0(3). (35) 

a=l 

0(3) stands for cubic or higher order terms. In one- loop approximation it is sufficient that one 
keeps only terms up to quadratic powers in the fluctuation fields which are identified here with 
the Vi and pa degrees of freedom. However we have kept the quartic term Y^a=i Pa equation 
()35|) for other purposes which will become clearer shortly. The action S'Ar[Aa] is then given by 

3 

Sm[K] = SN-i[Da] + vfn^jVj + Y}(^pI + Mpt] + 0(3). (36) 



a=l 



where clearly 

SN-i[Da] = -^Tr^.^[Da, D^f + !^e^,,Tr^_i[D„, D^]D, + /STr^.^Dl + MTr^_,{Dl)\ 

(37) 

The operators Qij are given explicitly by 

Qij = (^2p + (2M + N)Dl^ 6ij + (2M - N)DiDj - 2ia'^Fij. (38) 

From equations (jHUI) and (jHHj) it is quite clear that the (A^ — 1)— dimensional vectors Vi play 
exactly the role of (bosonic) quark fields moving in the background of a covariant U{1)^~^ 
gauge field Da- On the other hand the logarithmic potential ( equation (j^Tj) ) takes the form 

Af-l 

log(AA - Xb? = Yl log(^^ - + Y - /^s)'- (39) 

A<B A<B A=l 

By integrating out fj, v* and pa we obtain the effective action S'^„ jDj, da\ given by 

X j dp3e-'3pi-^^p|+EL-^og{rfA-P3)^_ (40) 

The effective action reads therefore 

Sf.AD^, dj,] = S^_,m - Y ^og{dA - dsf - log [ /■ rfpe-^^^-^^^+S"- i^^^-^-^)^ 

A<B L 

+ Tr^Tr^-i log (^{2(3 + (2M + N)DD5^, + (2M - N)DiDj - 2ia^Fi}j . 

(41) 
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Tr2 is the 2— dimensional trace associated with the remaining U{1) rotational symmetry of the 
two matrices Di and D2 ( since is treated differently -diagonalized-in this approach) and 
Trjv-i is the usual trace over the matrices; here Di and D2 are (A^ — l)x(A^ — 1) matrices. 

The sum of the first two terms in (PT|) is nothing but the action (P7j) with the replacement 
Aa — >Da, Xa — ^dA and Tr^ — ^Tr^-i- Thus 

+ Tt^Ttn-i log (^(2/5 + (2M + A^)D^)5i, + (2M - N)DiD^ - 2ia''Fij 

(42) 

This action has U{1)^~^ gauge symmetry and a U{1) rotational symmetry since the matrix D^^ 
is diagonal. In the partition function the gauge symmetry can be easily enlarged to U{N — 1) 
by rotating the diagonal matrix ( back ) to a general form C3 given by D^, = UC^U'^ where 
U is an (A^ — 1) x (iV — 1) unitary matrix. We will have the simultaneous rotations = UCiU^ . 
The action SN-i[Di,dA\ becomes given by (fT^ with the replacement Xa — >Ca and where the 
trace is normalized such that Trl = N — 1. We write the above result in the following suggestive 
form 



Sf_^[D,,dA] = ^^_i[A,c^A]-log 



log 



dpe 



-l3p^-Mp^+Y.lzl\og{dA-p? 



+ 



Tr^TrN-i log (^(2/? + (2M + N)Dl)S,j + (2M - N)DiDj - 2ia^F,,^ . 



(43) 



I is precisely the one-loop contribution to the classical action (fTSjl coming from integrating 
out from the model only one row and one column. In the large N limit we can treat as a 
continuous variable and thus we can simply obtain the full one-loop contribution to the classical 
action ()13p by integration over of the above result. 



3.2 



The S%- 



to-Matrix phase transition 



We are interested in particular in verifying the stability of the fuzzy sphere ground state (jl5|) 
under quantum fluctuations. We consider therefore the background Da = (y4>La where La are 
the generators of SU{2) in the irreducible representation ^^^-^ which are of size (A^ — 1) x (A^ — 1). 
(j) is the field associated with the fluctuations of the radius R = a. The classical potential from 
fjH7|) is given by 



v[4>] 



N -2 
N + 1 



2c2a^ 



1 -|- m 



1 

3' 



+ 



A^-2 

ivTI 



+ 



A^-2(A^-1)2- 
iVTI A^2 _ 1 



m 
~2 



(44) 
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where we have used the relations a = ay/N, j3 = —a^m? and M = We have also 

kept terms of order A^^ only which dominates in the large limit. Since = aipL^ and 

becomes in the large N limit 



{L^)ab = ^A^AB where = ^ — A we have = a(j)mA- Thus the last integral in (j84j) 



The constant c is independent of the field (p. We also compute in the large N limit 

2P + {2M + N)DI = a^(-2m^ + ^'~^f (c2 + mV) 

A/^ — 1 

= a' { —ct>' - - + m202 _ + o(^) j . (46) 



(2M - N)DiDj 



_ 2A^ . 

A^2 _ 1 

Ar2 



[m -C2) 



^01 1" N , 2m2 



2L,; 



VN^ - 2N 



(47) 



-2ia ^/c^ 



m-2N 

1 



.,5^|iV-l-- + 0(^) 



2L, 



VN^ - 2N 



(48) 



The leading quantum contribution of the effective potential from and (jlEI)-(|lHl) is 

given by 



6v; 



eS 

N-1 



-(2iV - 1) log + TrsTr log 



4 v^*^ 

(2A^ — 3) log0 + a constant independent of 



(49) 



As we will check in appendix A all higher order terms in equations (jini)-(0HI) will give vanishingly 
small quantum contributions in the large A^ limit. Thus is the one- loop correction of the 
effective potential coming from integrating out one column and one row from the theory. The 
full one-loop correction coming from integrating out all columns and rows is obtained by a 
simple integration over A^. We get 



V, 



eff 
N-1 



{N^-3N) log0. 



(50) 



By putting and together we get the effective potential 
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eff 



1 + 



m 
T 



+ A^^ iog0 + O(A^). 



(51) 



Let us point out here that this potential was derived elsewhere using a completely different 
(much simpler) argument involving gauge fixing the original action ()13|) and then computing 
the full one-loop effective action in the background field method. The argument in this article 
is however superior from the point of view that it ( manifestly) preserves gauge symmetry at all 
stages of the calculation since we are not fixing any gauge in the usual sense [22]. See also [23]. 

It is not difficult to check that the corresponding equation of motion of the potential (j^T|) 
admits two real solutions where we can identify the one with the least energy with the actual 
radius of the sphere. This however is only true up to a certain value of the coupling constant 
a where the quartic equation ceases to have any real solution and as a consequence the fuzzy 
sphere solution (fT5|l ceases to exist. In other words the potential V^g below the value of the 
coupling constant becomes unbounded and the fuzzy sphere collapses. The critical value can 
be easily computed and one finds 



8(1 + m2) 



l + \ 1 + 



32m2(l 
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(52) 



and 



at 



[1 + m 



2\ i4 



m 



Extrapolating to large masses we obtain the scaling behaviour 



a* 



[- 



+ V2- 



(53) 



(54) 



In other words the phase transition happens each time at a smaller value of the coupling 
constant d and thus the fuzzy sphere is more stable. This one-loop result is compared to the 
non-perturbative result coming from the Monte Carlo simulation of the model (|13|) with the 
constraint (fTTj) in figure 1. As one can immediately see there is an excellent agreement. In this 
sense the one-loop result for the U{1) is exact. Let us finally report that this phase transition 
was also observed in 4 dimensions on S| x S|. See the first reference of [15]. 



4 The small one-plaquette model limit on 8% 

We have found in the one-loop calculation as well as in numerical simulation that the 
presence of the normal field ()2H) is what causes the model to undergo the above first order 
phase transition from the fuzzy sphere to a matrix phase where the fuzzy sphere collapses 
under quantum fluctuation. At the level of perturbation theory of the gauge field Aa this shows 
up in the form of a compact UV-IR mixing phenomena which goes to the usual singular UV-IR 

13 
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Figure 1: The phase diagram of the S^— to-matrix phase transition. The fuzzy sphere phase 
is above the sohd hne while the matrix phase is below it. In this figure is the Monte Carlo 
measurement of the critical value . 

mixing on the NC plane in some appropriate planar limit of the sphere. In the large m limit 
we also have shown that these two ( possibly related ) effects disappear [22]. 

As it turns out there is some signature in Monte Carlo simulation of the model (|13|) with 
the constraint ()17p for the existence of another kind of phase transition which seems to be 
unrelated to the S^-to-matrix phase transition and which generically persists even in the large 
m limit. This latter phase transition resembles very much the third order one-plaquette phase 
transition in ordinary two dimensional gauge theory. In particular the agreement in the weak 
regime between the simulation and the theory ( which we will present now ) is excellent. Let us 
also say that this transition starts to appear when the critical value as m increases becomes 
less than the value and it becomes more pronounced as a^, decreases further away from this 
value. The new phase transition thus occurs at 



Our goal in this section is to give a detailed theoretical model which describes this transition. 
This construction is motivated by [24,26]. 

We start by making the observation that in the large m — ^oo limit we can set $ = as one 
can immediately see from the action (|^^ and the partition function (j211)- Indeed we have for 
m — >oo 



In other words the normal scalar field $ becomes infinitely heavy ( m is precisely its mass ) 
and thus decouples from the rest of the dynamics. Hence we can effectively impose the extra 



3.35 



(55) 





(56) 
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constraint = 0^02 on the field in this hmit m — >oo. In terms of Da = La + Aa = ^ 
this constraint reads 

The action ( if (fTTj) is also satisfied ) becomes 

S = SYM + Scs + So + ^\og(^^]. (58) 



2 V 

The Yang-Mills and Chern-Simons-like actions are given respectively by 

Scs = -TT^Tr [tabcFabDc + " ^2] = --^TreabcFabDc wlicn m — ^00. (59) 

In above 9^ = and 5*0 = S[Aa = 0] = — |a^C2 — |a^C2m^ . In the continuum large N — ^oo 
limit the constraint = o?C2 becomes the usual requirement that the normal component of 
the gauge field on the sphere is zero, viz $ = Ua-Aa = 0. Moreover the Chern-Simons-like 
action vanishes in this limit by this same condition ^ = Ua-Aa = because it will involve the 
integral of a 3— form over a 2— dimensional manifold. Hence S = Sym in the large N — >oo 
limit provided we also impose the condition $ = 0. In summary if we take the limit m — >oo 
first and then we take the continuum limit N — >oo we obtain a U{1) action on the ordinary 
sphere, viz 



S-So = — log j+ — J -{^CaAb - ^CbAa + eabAcY- (60) 

By construction the continuum gauge field ( which should be easily distinguished from its 
corresponding operator on the fuzzy sphere although we are using the same symbol Aa for both 
quantities ) is strictly tangent. See appendix B for more detail. 

If we study instead the action in the limit m — >oo first and then in the continuum 
limit N — i>oo then we will find the action 

= j '^{iCaAb-lCbAa + eabcAcf. (61) 

Again the continuum gauge field Aa is strictly tangent. The only difference with the previous 
case is the extra piece (5'^0q which we pulled out from the action in the process of writing it 
only in terms of a tangent gauge field. 

We now relate this action with the one-plaquette action. To this end we introduce the 
2Nx2N idempotent 

'y = j^{l2N + 2aaLa) , 7' = 1 (62) 
15 



where a a are the usual Pauh matrices. It has eigenvalues +1 and —1 with multiplicities + 1 
and — 1 respectively. We introduce the covariant derivative Da = La + Aa through a gauged 
idempotent 7d as follows 

1 

7d = 7 



7 = ^(1 + 2aaD,) = 7 + ^a^A, , 7^ = 1 + ^$ + (63) 

Since we are interested in the large m — >oo limit we may as well set $ = in above. Clearly 
7£) has the same spectrum as 7. Thus it is a continuous deformation of 7 in the sense that 
there exists a U{2N) unitary transformation U such that 7/) = U'~fU~^. Furthermore if U — >UT 
where Te.U{N + 1) or TeU{N — 1) then 7 — >T'yT^ = 7 and as a consequence 7d — ^7d- So 
71) is an element of the (iAr— Grassmannian manifold U{2N)/U{N + l)xU{N — 1). We compute 
the dimension c^tv as follows 

dj^ = AN^ -{N+ 1)2 - (A^ - 1)2 = 2N^ - 2. (64) 

This is exactly the correct number of degrees of freedom in a gauge theory on the sphere without 
normal scalar field or with a normal scalar field frozen to some fixed value. The 2 counts the 
zero modes which decouple because of commutators. 

The original U{N) gauge symmetry acts on the covariant derivatives Da as 

Dl = gDag^ , geU{N). (65) 

This symmetry will be enlarged to the following U {2N) symmetry. We introduce a tentative link 
variable W ( a 2Nx2N unitary matrix )hj W = 77^- The extended U{2N) symmetry will then 
act on W as follows W — >VWV^ , V&U{2N). It is clear that this transformation property of 
W can only be obtained if we impose the following transformation properties 7 — >V'~fV^ and 
7£) — >V'-f£)V~^ on 7 and 7^ respectively. Hence the U{N) subgroup of this U{2N) symmetry 
which will act on Da as Da — >gDag~^ will also have to act on La as La — >gLag~^, i.e 

= gDag"- , LI = gLag^ , geU{N). (66) 

It is not difficult to see that the two sets of gauge transformations ()65p and ()66p are identical if 
we are looking at the action (fT^ since it only depends on Xa = aDa and not on La. However 
for the gauge field Aa there is certainly a difference between the two sets of transformations 
and (jnni)- Under (jH^ we have Aa — >gAag^+g[La, g~^] whereas under (jHUj) we have Aa — >gAag^ ■ 
The actions as written in and (j^Hj) are invariant only in the first case. 

We want thus to modify the definition of the link variable W so that we have ()65p and not 
fj66p . In other words under this new definition the fixed background La will not rotate whereas 
the gauge field Aa will transform correctly as Aa — >gAag^ + g[La,g^]- Towards this end we 
introduce another covariant derivative D'^ = La + A'^ through the gauged idempotent 7^' given 
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by 

Id' = 1 

1 2 8 /c~ 2 

7 = ^(l + 2aX) =7 + ^^a< , = 1 + ^^' + j^eab,aXb- (67) 

As before we will also set $' = 0. From the two idempotents 7^) and 7^' we construct the link 
variable W as follows 

W = -i^nD. (68) 

The extended U{2N) symmetry will then act on W as follows 

W — >VWV+ , VeU{2N). (69) 

This transformation property of W can only be obtained if we impose the following transfor- 
mation properties 7^,' — >V'j^'V~^ and — ^V'JdV'^ on 7^' and •jd respectively. Hence the 
U{N) subgroup of this U{2N) symmetry which will act on Da as Da — >gDag~^ will also have 
to act on D'^ as D'a — >gD'^g^, i.e 

= gDag^ , = gD[g^ , geU{N). (70) 

Under these transformations the gauge fields Aa and A'^ transform as Aa — >gAag^ + g[La, g'^] 
and — >gX^g^ + glLa^g^] respectively like we want. 

Remark also that for every fixed configuration X^ the link variable W contains the same 
degrees of freedom contained in 7^3. To see this we will go to the basis in which 7^' is diagonal, 
viz 

- ^ ( '7' ) ■ 

In this basis '')d and W will have the following generic forms 

Wi = is an (A^ + 1) x (A^ + 1) matrix, W2 = is an (A^ - 1) x (A^ - 1) matrix and Wu 
is an (A^ + 1) x (A^ — 1) matrix whereas the hermitian adjoint is an (A^ — 1) x (A^ + 1) 
matrix. Since 7I) = 1 or equivalently W^W = 1 we must also have the conditions 

W^Wi + Wi2Wt2 = 1 

W^W2 + W+Wu = 1 

WiWu + W12W2 = 0. (73) 

Knowing W12 will determine completely the matrix W (or equivalently 'Jd ) and hence we have 
2(A^ + 1)(A^ — 1) = 2A^^ — 2 degrees of freedom which agrees with 
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4.1 The coordinate transformation {Ai, A2, A^) — >{W,^) 

The main idea is that we want to reparametrize the gauge field on S| in terms of the fuzzy 
hnk variable W and the normal scalar field $. In other words we want to replace the triplet 
{Ai, A2, A3) with {W, $). It is the link variable W which contains the degrees of freedom of the 
gauge field which are tangent to the sphere as is shown by the result (jMj) . Thus in summary 
we have the coordinate transformation 

{A,,A2,As)-^{W,<l>) (74) 

First we need to show that we have indeed the correct measure. Namely one must show that 
we have 

j dAidA2dA^ j ^^^^ 

where cat is some constant of proportionality which can only depend on A^. In order to compute 
the measure we will compute the quantity Tr2N{dW)^ dW where Tr2N denotes the 2N x 2N 
dimensional trace. For this exercise the scalar field $ will not be assumed to be fixed whereas 
the other gauge configuration and its corresponding normal scalar field are supposed to 
be some constant backgrounds. From the definition W = 7£)'7_d and equations (jU^ and ^7\i 
one can easily compute 



or equivalently 



2 2 2 2 1 

W = 1 + —-faaAa + —aaA'-f $ $' + 0(— ) (76) 



dW = ^^aadAa - ^d^ + O(^). (77) 



Hence a straightforward calculation yields the measure 

Tr2N{dWydW = ^Tr{dA,f + ^Tr{d^f - ^^Trd^diL^A, + A,L,) + O(^). (78) 
By using the identity 2^/c2$ = AaLa + LaAa + we arrive at the result 

Tr2N{dWYdW = -^JridA^Y - ^Tr{d<l>r + O(^) (79) 

The correct ( more suggestive ) way of writing this equation is the following 

1 

Tr{dAaf = —Tr2N{dW)+dW + Tr{d^f + O(-). (80) 
8 iV 

In the large — >oo limit it is obvious that this equation implies (j75j) which is what we desire. 
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4.2 The U{1) gauge action as a linear one-plaquette model 

It remains now to show that the enlarged U{2N) symmetry reduces to its U (N) subgroup 
in the large N limit. The starting point is the 2 A^— dimensional one-plaquette action with a 
positive coupling constant A, viz 

Sp = —Tr2N{W + W+ -2) (81) 
A 



with the constraints 



We have the path integral 



W = 7z,'7D , $ = ^^It^- (82) 



c/7^'c/$'5($') f dWd^5{^)e^p. (83) 



c% is the constant which appears in (ffH|l . The extra integrations over 7^5/ and $' ( in other 
words over D'^ ) is included in order to maintain gauge invariance of the path integral. The 
integration over W is done along the orbit W = 7£)'7d inside the full U{2N) gauge group. In 
above we have also to integrate over configurations Da and such that $ = and $' = since 
we are only interested in the limit m — >oo of the model ()23j) . Furthermore we can conclude 
from the result (|7H)l that in the large N limit this path integral can be written as 

dX^5{^') j dAa5{^)e^''. (84) 

We need now to check what happens to the action Sp in the large N limit. This is done in 
appendix C and one finds 



N 



(85) 



The constraints = C2 and D'^ = C2 ( or equivalently $ = and $' = ) become in terms of 
the variables Aa = \Aa — \ ^'^d Sa = Da + D'^ ( or equivalently Aa = ^Aa + ^A'^ ) 

Si + AAl = Ac2 ^Al + {La, Aa] + Al = ^ 

[Sa, Aa} = Q^ {La, Aa} + {ia, ^} = 0. (86) 

In the continuum limit these two constraints becomes UaAa = and riaAa = respectively. By 
using the first constraint we can rewrite the action in the form 



S, 



N 
J 



— Tri^L, + ^Tri^ + — Tr U[La, A,] - t[L,, Aa] + e,,,A J + 0{ 



1 , 



^7) 
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The leading contribution in the action Sp as written in equation ()85p is a simple Gaussian which 
is clearly dominated by the configuration Aa = 0. As a consequence the full path integration 
over Aa is dominated hj Aa = 0. This yields a zero action which is obviously not what we want. 
Furthermore the path integration over Aa diverges since this action (jH3j) does not depend on 
these matrices. On the other hand the Gaussian term becomes in equation ()87|) ( after using 
the constraint ) a quadratic integral over the matrices Aa but with a wrong sign since the first 
term converges to in the limit ( see appendix B ). Thus the path integration over the three 
matrices Aa will again diverge. The one-plaquette action Sp by itself is therefore not enough 
to obtain a U{1) action on the sphere in the continuum large limit. 



4.3 A quadratic one-plaquette action 

Towards the end of constructing a U{1) action on the fuzzy sphere using the one-plaquette 
variable W we add to the action Sp the following quadratic one-plaquette action ( where A' is 
the corresponding coupling constant ) 



S' 



^^TrMW + W^' - 2) 
A 



Remark the extra minus sign in front of this action, i.e A is a positive coupling constant. As 
before we need now to compute the large A^ limit of this quadratic one-plaquette action. This 
is also done in appendix C and one finds the result 



s: 



N 



'-^TrAaLa + ^TrAl + ^Tr(i^ + {i., L J)^ + ^Jr [i[La, A, 



- i[Lb,Aa]+eabcAc] +0 



(89) 



By putting the one-plaquette actions (fHTjl and (jHHI) together we obtain the total one- 
plaquette action 



Sp + Sp = -^Tr{{Aa,La} + Al)-^Tr{{Aa,La} + AlY-^Tr{z[La,A,] 



t[L,,Aa]+eatcA,] +0{^,)-0{^). 



(90) 



The positive coupling constant Ai is defined in terms of A abd A' by 

1 _ 1 4 



(91) 



The effect of the dominant terms ( the first two terms in the above action ) is now precisely 
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what we want. The path integral over the three matrices Aa is given by 



= 2^""" j dAa5(^^{Xa,Aa} + ^{Aa,Aa}y(^-^^ exp 



In the large N limit the first term in the exponent dominates ( see below ) and as a consequence 
the path integral over the three matrices Aa becomes a simple Gaussian. Since the second 
constraint inside the integral has the effect of reducing the number of independent matrices Aa 
to just two we obtain the final result 

S(A,) ^ 2^'^s(^lx,,A,}\ JdAj(^^lx,,A,}\ exp[-^TrAl\ 

Another ( more correct ) way of understanding this result is to note that this path integral is 
dominated in the large N limit by the configurations Aa = 0. 

The path integral of the one-plaquette model ( with an action Sp + S'p ) becomes in the 
large N limit as follows 

Z'p= [ dAaS(l{xa,Aa}]e'''^ (94) 



2 



where 



Sf = NHog{^)-^^Tr(^^[La,A,]~^[U,Aa]+ea,.A)j\o{^^^ 

(95) 

Notice that this action is invariant not only under the trivial original gauge transformation 
law Aa — >Aa but also it is invariant under the non-trivial gauge transformation Aa — >Aa + 
g[La,g'^] where gE.U{N). This emergent new gauge transformation of Aa is identical to the 
transformation property of a U{1) gauge field on the sphere. Therefore the action 5*1? given 
by the above equation is essentially the same U{1) action —{S — Sq) given in equation (plj) 
provided we make the following identification 

16 _ 16 1 A I _ 
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between the U{1) gauge coupling constant g on the fuzzy sphere and the one-plaquette model 
coupling constant Ai. The action becomes 

(97) 

Let us remark that in this large limit in which g is kept fixed the one-plaquette coupling con- 
stant Ai goes to zero. Hence the fuzzy sphere action with fixed coupling constant g corresponds 
in this particular limit to the one-plaquette gauge field in the weak regime and agreement be- 
tween the two is expected only for weak couplings ( large values of a ) . To see this more clearly 
we notice that in terms of A and A' the limit Ai = , — >0 is equivalent to the limit A — >0 
for fixed A or vice versa, i.e to the limit A — >0 for fixed A. Furthermore Ai going to is also 
equivalent to the limit when both A and A' go to zero. Clearly all these possibilities correspond 
to the one-plaquette gauge field in the weak regime. 

Finally we remark that the constant term in ()6U|) depends on the mass parameter m. Thus 
by comparing between the constant terms in (jUUI) and (j^7|) we can determine m? as a function 
of g^ ( or equivalents ) and N. We find = = ^^0^. 

4.4 The one-plaquette path integral 

Instead of (j83|) we will therefore consider in the remainder of this article the following ( 
corrected or generalized ) one-plaquette path integral 

Zp = c% I d-fjj'd^'S{^') I dWd^5{^)e^p^^'p. (98) 

In analogy with (f?^ we decompose the 2Nx2N matrices 7£)'7z) and W as follows 

In particular Wi = is an (A^ + 1) x (A^ + 1) matrix, W2 = W2 is an (A^ — 1) x (A^ — 1) matrix 
and W12 is an (A^ + 1) x (A^ — 1) matrix whereas the hermitian adjoint is an (A^ — 1) x (A^ + 1) 
matrix. Since W^W = 1 we have the conditions 

W^Wi + Wi2W^2 = 1 . W+W2 + W+Wu = 1 , WiWu + W12W2 = 0. (100) 
First we observe that in this basis the metric becomes 

Tr2N{dW)+dW = TrN+iidWi^dWi + TrN-i{dW2)^dW2 

+ TrN+idWuidWu^ + TrN_i{dWuVdWi2. (101) 
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Hence we can immediately conclude that the path integral over W can be rewritten ( by 
neglecting an overall proportionality factor ) as 



I dWoc I dWi2 j dWt2 f dWi j dW2. 



(102) 



Furthermore we can show that in this basis the actions Sp and Sp take the form 



Sp = —TrN+^{Wi + W+-2) + —TrN,^{-W2-W+~2) (103) 
A A 



and 



S'p= - yTrN+i{W^ + W^^-2)-yTr^^,{W^ + W+^-2) 

2N 2N 
+ —TrN+,WuW+ + —TrN_,W+W,2. (104) 

Thus the off-diagonal matrices W12 and ( as opposed to the diagonal matrices Wi and W2 
) appear only in the action S'p. 

Let us recall that since the integration over W is done along the orbit W = 7£)'7d inside 
U{2N) and since in the large limit both 7^' and 7/3 approach the usual chirality operator 
7 = UaCTa we sce that W approaches the identity matrix in this limit. It is in this sense that W 
yields in the continuum large limit a small one-plaquette model. 

Let us now explain how we will approximate the above path integral in the continuum 
large A^ limit. From one hand we have the following limiting constraint W = 7£i'7d — ^'^2N^ 
which means that when A^ — ^00 we have the behaviour Wi = (7£)'7d)i — ^Ia^+Ij —W2 = 
— (7^/7^)2 — ^Iat-i and W12 = (7£)'7d)i2 — ^0. From the other hand since W must be always a 
unitary matrix and since the off-diagonal parts Wu and tend to zero the matrices Wi and 
—W2 become in this approximations (A^ + 1) x (A^ + 1) and (A^ — 1) x {N — 1) unitary matrices 
respectively ( which are close to the identity ) in accordance with equations ()100|) . Let us also 
stress the fact that the strict limits of Wi, —W2 and W12 are independent of 7^'. For example 
the matrix Wi goes always to the same limit 1 at+i for all matrices 7^)' . 

The main approximation which we will adopt in this article consists therefore in replacing 
the constraint W = 7£i'7d with the simpler constraint W — ^l2N by taking the diagonal parts 
Wi and —W2 to be two arbitrary, i.e independent of 7^', unitary matrices which are very close 
to the identities Iat+i and Iat-i respectively while allowing the off-diagonal parts Wu and 
to go to zero. We observe that by including only Wi and —W2 in this approximation we are 
including in the limit precisely the correct number of degrees of freedom tangent to the sphere, 
viz 2A^^. Thus in this approximation the integrations over $, $' and 7^' decouple while the 
integrations over W12 and are dominated by W12 = = 0. There remains the two 
independent path integrals over Wi and —W2 which are clearly equal in the strict limit since 

^Notice that W = jj^'-fo goes to I2N independently of any basis. 
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the matrix dimension of Wi approaches the matrix dimension of —W2 for large A^. Thus the 
path integral Z'p reduces ( by neglecting also an overall proportionality factor ) to 



Zp cx [Zp{X,X')]^ (105) 

where 

Zp{\, a') = j dWi exp |yTr(iyi + - 2) - yTr(iy2 + |y+2 _ 2)|. (106) 

The path integral of a 2— dimensional U{N) gauge theory in the axial gauge Ai = on a 
lattice with volume V and lattice spacing a is given by Zp{\, 00)^/"' where Zp{X, 00) is the 
above partition function ()106p for A' = 00, i.e the partition function of the one-plaquette model 
Sp = ^Tr{Wi + W^ — 2). Next we need to understand the effect of the addition of the quadratic 
one-plaquette action 5*^ = —YTr{Wi+W^'^ — 2). Formally the partition function Zp{\, X')^^"''^ 
for any value of the coupling constant A can be obtained by expanding the model Sp + Sp 
around A' = 00. Thus it is not difficult to observe that the one-plaquette action Sp + S'p does 
also lead to ( a more complicated ) U{N) gauge theory in two dimensions. The U{N) gauge 
coupling constant gf is simply given by 

\ = Na'\. (107) 
9i 

Therefore we can see that the partition function Z'p of a U{1) gauge field on the fuzzy sphere is 
proportional to the partition function of a g'enerafeec? 2— dimensional U{N) gauge theory in the 
axial gauge = on a lattice with two plaquettes. This doubling of plaquettes is reminiscent 
of the usual doubling of points in Connes standard model. The U{1) gauge coupling constant 
g"^ and the U{N) gauge coupling constant gf are related ( from (jHS)) and ()107|) ) by the equation 

'--1 . (108) 



g\ 64 \X 

It is quite natural to require the two coupling constants g^ and g\ to be equal which means we 
must choose the lattice spacing a such that 

4 64 A' , , 

4.5 Saddle point solution 

We are therefore interested in the A^— dimensional one-plaquette model 

Zp{\,\') = j dWexp(^^Tr{W + -2) - ^Tr{W^ + W+^ -2)y (110) 
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Let us recall that dW is the U{N) Haar measure. We can immediately diagonalize the link 
variable W by writing W = TDT^ where T is some U{N) matrix and D is diagonal with 
elements equal to the eigenvalues exp{i6i) of W. In other words Dij = 5ijexp{i0i). The 
integration over T can be done trivially and one ends up with the path integral 

^p(A,A') = j X^^^de.e'''^. (Ill) 

The action Sn contains besides the Wilson actions ^Tr{W + — 2) = ^ Yl!i=i cos9i — ^ and 
-^Tr(VF^ + W'^'^ — 2) = Y cos29i — contributions coming from the usual Vandermonde 
determinant. Explicitly the total action reads 



= ^5;cos«.-f^cos2«, + i^^ln(=inM^?^^ + f. (112) 

In the large limit we can resort to the method of steepest descent to evaluate the path 
integral Zp(A, A ) . The partition function will be dominated by the solution of the equation 
= which is a minimum of the action S]\j. Before we proceed to the solution we need to 
take into account the following crucial property. Since the link variable W tends to one in the 
large — *-oo limit we can conclude that all the angles 6i tend to in this limit and thus we 
can consider instead of the full one-plaquette model action ()112|) a small one-plaquette model 
action by including corrections up to the quadratic order in the angles 9i. We obtain 



^2 . — • 



A2 is given by 



111 , , 



For the consistency of the solution below the coupling constant Ai must be negative ( as opposed 
to the classical model where Ai was assumed positive ) and as a consequence the coupling 
constant A2 is always positive. As it turns out most of the classical arguments of sections 4.2 
and 4.3 will go through unchanged when Ai is taken negative. 

Thus in the following quantum theory of the model we will identify the effective one- 
plaquette action S'ff with the fuzzy sphere action S — Sq { which is to be compared with 
the classical identification —S'j^ = S* — Sq) and hence we must make the following identification 
of the coupling constants 

' (115) 



N'^Xi 4^2 4JV2 
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This is precisely due as we have said to the fact that Ai becomes negative in the quantum 
theory. In the continuum large N limit where is kept fixed instead of Ai we can see that ^ 
scales with A^^ and as a consequence 

The saddle point solution must satisfy the equation of motion 



I sin Oi-^ sin = 1 y cot (117) 
A A' 2 ^ ^ 



The equation of motion ()133p takes ( in the limit — >-oo when all the angles tend to zero 
) the form 



In order to solve the above problem we introduce the potential V{{6i}) defined through its first 



6i, i = 1, ...,N. The trace uj{z) of the resolvent of M is given by 



derivative — ^Jg ''' =V {d.i) = ^ and also the N x N matrix M defined through its eigenvalues 



, , 1 1 1 1 , , 

i 

The condition ()118|) can then be rewritten as follows 

.^z) - ^.{z) + V'izMz) = -R{z) ^-^Y: ^'^'IZl^^'^ - (120) 

i 

In the large limit we can also introduce a density of eigenvalues p{6) which is positive definite 
and normalized to one ; p{9) > 0, J d6p{0) = 1 [ N p{6) is the number of eigenvalues in the 
range [6 — dO /2,6 + dO /2] ]. Thus the sum will be replaced by Yli = ^ / d9p{9) and one obtain 

uj\z) + V'{z)uj{z) = -Riz) = - r d9p(9) ^ (^) W _ (121) 

The trace of the resolvant is now given by 

u{z) = r d9pi9)-^. (122) 

J-e, 0- z 

The density of eigenvalues p{9) should satisfy d9p{9) = 1 and p{9)>0 for all angles 
—9^<9<9^:. We can easily solve this problem since we can compute 

R{z) = ^ (123) 
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and 



a{z) = ^V2A2-^2. (124) 



The solution of the equation of motion is immediately given by 

uj^{z) = -\v\z)±\a{z) 

z 



±±^2\2-z\ (125) 



The function 00 {z) is a multi- valued function of z with branch points at z = ±zo = ±a/2A2. 
Since the potential V has only one minimum at = the density of eigenvalues must have only 
one support centered around this minimum. This support is clearly in the range between —zq 
and +zq. In terms of the resolvent u^z) the density of eigenvalues is defined by 

piz) = ^(^ + ^^)-^(^-^^). (126) 
2-Ki 

Lj{z + ie) is the trace of the resolvent of M computed with a contour in the upper half complex 
plane and we choose for it the plus sign, viz u^z + ie) = uj+{z). Similarly u^z — ie) is the trace 
of the resolvent of M computed with a contour in the lower half complex plane and we choose 
for it the minus sign, viz uj{z — ie) = u!-{z). We obtain therefore 

p{0) = -^V2X^. (127) 

It is obvious that this density of eigenvalues is only defined for angles 6 which are in the range 
—\/2X2<0<\/2X2. However the value of the critical angle 6^, should be determined from the 
normalization condition J^* d6p(6) = 1. This condition yields the value 

= ^/2\2. (128) 

4.6 The one-plaquette phase transition 

It is quite obvious that the action pi3|) is an excellent approximation of ()112j) for all angles 
9i in the range 

- \<(^^<\- (129) 

The particular value | comes from the fact that the expansion of the quadratic one-plaquette 
action S'p will converge to the original expression only for small 9i in the above range. The 
expansions of the linear one-plaquette action Sp and of the Vandermonde action will converge 
to the original expressions for 9i in the range — l<^i<l. 
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The solution p27|) with the critical angle p28p is then valid only for very small values of 
the coupling constant A2. Indeed it is only in this regime of small A2 where the fuzzy sphere 
action with fixed coupling constant g is expected to correspond to the one-plaquette model as 
we have discussed previously. However in order to find the critical value of A2 we need to extend 
the solution ()127|) to higher values of A2. To this end we note that the action ()113|) can also be 
obtained from the effective one-plaquette model 

Sf = ^Tr{W,s + - 2) 
^2 

For small in the range 

- l<^f <1 (131) 
The total effective one-plaquette action becomes 

= -^j:i(^fr+^j:^- ^ +o((rY)- (132) 

The action ()132|) must be identical to the action ()113|1 and hence we must have {Oi^Y = -x^^^f- 
From the two ranges (HSni) and (USD) we conclude that Of = 281 and \ f = 4A2. 
The saddle point solution of the action (|13U|) must satisfy the equation of motion 

<-) -1 ncR neS 

J_ Sin 9f = - y cot ^ (133) 

In the continuum large N limit this equation becomes 



^sin6'efr= / (ircffp(rcff) cot ^. (134) 



A. 

By using the expansion cot = 2 Xl-^i ( sin cos nr — cos n9 sin nr) we can solve this equa- 
tion quite easily in the strong-coupling phase ( large values of A2 ) and one finds the solution 



p{Ocs) = 77- + — TrfT cos6'efr. (135) 
27r 71X2 



However it is obvious that this solution makes sense only where the density of eigenvalues is 
positive definite, i.e for A^''^ such that 

^ - ^>0 ^(Af )* = 2 ^A; = 0.5. (136) 
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This strong-coupling solution should certainly work for large enough values of However this 
is not the regime we want. To find the solution for small values of A2 the only difference with 
the above analysis is that the range of the eigenvalues is now [—6^, +6^] instead of [— vr, +7r] 
where 9^ is an angle less than vr which is a function of A2. It is only in this regime of small A2 
where the fuzzy sphere action with fixed coupling constant g is expected to correspond to the 
one-plaquette. In the strong regime deviations become significant near the sphere-to-matrix 
transition. Finding the solution in the weak-coupling phase for the effective action (jl32p is a 
more involved exercise. This is done in [25] with the result 

p(M = ^cos-y— -sm -. (137) 



sm- = y— . (138) 

It is very easy to verify that the this density of eigenvalues and critical angle will reduce to the 
solution (jl27p and the critical angle (jl28j) when the angles are taken to be very small. 

The above computed critical value A2* = 0.5 leads to the critical value of the coupling 
constant a 

64 

= — = 128 ^ a, = 3.36 (139) 

which is to be compared with the observed value 

a, = 3.35±0.25 (140) 

Indeed for U{1) theory we observe in Monte Carlo simulation of the model (fT^ with the relation 
(fT7|) the value = 3.35±0.25 ( or equivalently the value A2* = 0.51±0.15 ). Indeed for very 
large values of the mass parameter m we observe two critical lines ( see figure 2); the lower line 
is the S^-to-matrix critical line discussed previously. This line comes from the measurement 
of the critical value as = a* from the action. The upper line is the one-plaquette critical line 
which we can fit to the curve 

a„ = \^r\ ^ + 3.35±0.25. 141 

L ^2 J 

Remark that this curve saturates in the limit m — >oo around the value 3.35. The points 
on figure 2 comes from the measurement of the position a* = ^J~N a^, of the peak in the specific 
heat which for large values of the mass captures the one-plaquette phase transition. For even 
larger values of m the peak disappears and in this case Ofp measures the position where the 
specific heat jumps discontinously to the value 1. 
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Figure 2: The phase diagram of the one-plaquette phase transition. 



4.7 The specific heat and effective potential in expansion 

We are now in a position to compute the quadratic average Q defined by the equation 



N 

i 

We obtain in the fuzzy one-plaquette solution ()127|1 the result 



Q = ^E^^ (142) 



Q = I d9p{9)9' = ^ (143) 



We need also to compute the non-local average 



2 



iln^ + ^. (144) 
2 2 2 ^ ' 



Si is a constant of integration given explicitly by 5*1 = ^ J^-^ dx\/\ — J^-^ dy\/l — ip- In ( 
. The action (|113|) is therefore given by 



-l + iln^ + ^. (145) 
2 2 2 2 ^ ' 
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Let us recall from equations ()105|) and p06p that we have actually two identical one-plaquette 
models and hence the above action must be multiplied by a factor of 2. Furthermore by 
comparing between and we can see that S'fF must be identified with NS^ ( or twice 
as much due to the above factor of 2 ) whereas we have found that the action 5* — 5*0 on the 
fuzzy sphere must be identified in the quantum theory with 5*1?. In other words the effective 
action on the fuzzy sphere is given by 

S = So + N^{-l + \n^ + S,) 

1 1 N'^ 

= --a%--a%m^-N'\na'' + N\-l-\n— + S,). (146) 
2 62 

In above we have also used equation ()116|1 . It is interesting to compare this effective action with 
the original effective action (j51|) obtained in the one-loop. If we set = 1 in (j51|) then we will 
find the same classical action as in the above equation, namely —\a^C2 — \ci^C2m?. However 
the quantum correction in (j3T|) in terms of a is by inspection given by A^^ In a which is different 
from the quantum correction in the above equation which is equal to — 4iV^ln(5. We also note 
that in the large m — >oo, then large — >oo limit the above action will be dominated by the 
classical mass-dependent term —}^a.'^C2m? ■ This is precisely what we observe in Monte Carlo 
simulation. See figure 3. 



The action for different values of m for N=8 
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Figure 3: The action for non-zero mass. The fit is given by the second term of equation ((T 



Finally we need to compute the specific heat. Towards this end we implement the scaling 

s 

T 



transformations S — and — The specific heat is then defined by 
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A straightforward calculation yields the very simple result 



:i48) 



Again this is what we observe in our numerical simulation of the U{1) gauge field on the fuzzy 
sphere in the weak regime. In the strong regime deviations are significant near the sphere-to- 
matrix transition. See figure 4. In this regime of strong couplings the action and specific heat 
are computed using the distribution of eigenvalues (jl37|) . We find 



S 



-L _4 J- _4 2 

--a C2 a C2m 

6 2 



1 

2 ' 



a , 
128' 



64 ^ 



and 



. a , 
128' 



(149) 



(150) 



We observe then that in the weak regime the specific heat is essentially given by 4c2 = A^^ — 1 
within statistical errors whereas in the strong regime the data does only follow the theoretical 
one-plaquette prediction away from the S^-to-matrix transition. This is presumably due to 
the effects of the matrix phase which becomes strong near the critical S^-to-matrix transition. 
Remark that the minimum of the specific heat is where the S^-to-matrix transition happens 
for large values of m. 



m =100,70 for N=12,16 



ra =70,N.12 
m''s1 00,N=16 
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Figure 4: The specific heat for very large values of m. 



Let us comment further on the quantum effective potential for U{1) gauge field on the 
fuzzy sphere S% in this 1/A^ expansion of the one-plaquette model. As we have said before by 
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comparing the effective potential ()146|) with the one-loop effective potential ()5H) in which 
we set = 1 we can see that the classical contribution is the same in both potentials whereas 
the quantum correction in (|5ip in terms of a is given by A^^ In a which is different from the 
quantum correction — 4iV^ In a in ()146p . This observation allows us to rewrite ( or to guess that 
) equation ()146|) ( should be rewritten ) in terms of the radius as follows 



S 



1 + m 



1 

3' 



m 

T 



-4A^2 log. 



(151) 



Now in contrast to what we have done so far in this article we will choose the mass parameter 
m to be proportional to A^. The simplest most natural choice is m = ffiN . The effects of the 
large mass limit will then be included implicitly in the continuum large N limit. This is in fact 
what was done in [24]. The above effective potential becomes 



-4Ar2 iQg, 



(152) 



It is very easy to verify that this potential admits a local minimum for all values of the coupling 
constant a. The minimum value (p^i^ is found to be given by 



mm 



1 + + 
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(153) 



In the limit m — >oo first and then — >oo ( considered in this article ) we can see that fh — >oo 
and hence 0min = 1 as expected. The fuzzy sphere ground state (fT^ is extremely stable in this 
limit and the the 1/A^ potential S ( as opposed to the one- loop potential ) is completely 
insensitive to the S^— to- matrix phase transition. 



5 Conclusion 

In this article we have shown explicitly that quantum noncommutative U{1) gauge field on 
the fuzzy sphere is equivalent ( at least in the fuzzy sphere-weak coupling phase ) to a 
quantum commutative 2— dimensional U{N) gauge field on a lattice with two plaquettes. 

By using the structure of the fuzzy sphere we have constructed a 2Nx2N matrix W given 
by equation (68) which was shown to contain the correct number of degrees of freedom tangent 
to the fuzzy sphere, namely 2A^^ degrees of freddom. The other N"^ degrees of freedom are 
contained obviously in the normal scalar field Indeed we have shown that the gauge field 
[Ai, A2, A^) is equivalent to (W, $) and dAidA2dA3 oc dWd^. The fuzzy sphere action (22) ( 
or equivalently (13) ) written in terms of Aa can therefore be rewritten in trems of W and ^. 
We have shown explicitly that in the limit where we can set $ = this action (22) and the 
action Sp + S'p given by (81) + (88) will tend to the same continuum limit, viz a U{1) gauge 
theory on the ordinary sphere. As a consequence the partition function of the fuzzy U{1) model 
in the limit m — >^oo can be given either by the large m limit of (24) or by equation (98). 
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Indeed the fuzzy partition function Z'p given by equation (98) is our starting point. It is 
found to be proportional to the partition function of a U{N) model in the axial gauge Ai = 
on a lattice with two plaquettes given by equations (|1(J5|) and (jl(J6|) . We remark that the U{N) 
theory consists of the canonical one-plaquette Wilson action Sp = ^Tr{Wi + Wi — 2) plus a 
novel quadratic one-plaquette action S'^ = — ^Tr(iyf + W^"^ — 2) together with the canonical 
measure dW . This is in fact the reason why W is called a link variable. The quadratic term was 
needed in order that the fuzzy sphere one-plaquette path integral Z'p converges to the sphere 
path integral in the large N limit. Remark also that the effective actions Sp and S'^ involve the 
NxN link variable Wi as opposed to the original actions Sp and S'p which involve the 2Nx2N 
link variable W. 

The doubling of plaquettes is a natural consequence of the model and it is reminiscent of 
the usual doubling of points in Connes standard model. However the ( other kind of ) doubling 
of U{1) fuzzy gauge fields was needed in order to have a gauge invariant formulation of the 
fuzzy one-plaquette model. In fact a covariant plaquette variable W can only be constructed 
out of two such U{1) fuzzy fields. 

The main results are given by equations (jl(J5|) and (jl(J6j) . It is therefore of paramount 
importance to find a more rigorous derivation of these two equations. Furthermore it will be 
very interesting to show that the large m limit of the path integral (24) and the path integral 
(98) are also equivalent for finite A^. Their large N equivalence used in this article is confirmed 
in our Monte Carlo simulation of the model ( which uses (24) ) where the measurement of the 
critical line a^, = 3.36 and the specific heat can be understood in very simple terms using the 
limit (105) of (98). In particular the value a^, = 3.36 seen in the simulation is precisely the 
Gross- Wadia-Witten one-plaquette 3rd order transition point as calculated in this article from 
the path integral (105). 

Since the plaquette variable W is small, i.e it approaches I2N in the large limit, we were 
able to show that the model in this limit reduces to a simple matrix model and as a consequence 
was easily solved. We computed the critical point and showed that it agrees with the observed 
value. We computed also the quantum effective potential and the specific heat for f/(l) gauge 
field on the fuzzy sphere S"^ in the 1/A^ expansion using this one-plaquette model. In particular 
the specific heat was found to be equal to 1 in the fuzzy sphere-weak coupling phase of the 
gauge field which agrees with the observed value 1 seen in Monte Carlo simulation. The value 1 
comes precisely because we have two plaquettes which approximate the noncommutative U{1) 
gauge field on the fuzzy sphere. In the fuzzy sphere-strong coupling phase deviations were 
found to be significant near the S^— to- matrix critical point. It will be very interesting to be 
able to extend this one-plaquette model to these large values of the gauge coupling constant g^, 
i.e to small values of a. The key to this we believe lies in improving the basic approximations 
of this article given in equations ()105|) and ()106|) . 

The most natural generalization of this work should include fermions in two dimensions [30] 
and as a consequence take into account topological excitations [29]. The best example which 
comes to mind is the Schwinger model [31]. Then one should seriously contemplate going to 4 
dimensions with the full might of QCD. Early steps towards this larger goal were taken in the 
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first reference of [15]. First we need to have a complete control over the phase diagram of the 
pure gauge model considered in this article [12]. 
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A Next-to-leading correction of the effective potential 

The next-to-leading contribution ( coming from the terms of order N and of order 1 in 
equations (PBj) - (PHj) ) is given by 



Tr2TrN-i log 



2 4m\., X 1,^ Ai . 1,, 1, 



(154) 



The 0— dependence is only in the second and third terms inside the logarithm. Let us also 
remark that the inverse of the operator 

Lij = aSij + bxiXj (155) 

is given by 

{L = —Sij H T^j- (156) 

a a X3-I-I 

lia = —b=l — j^ — then we can see that the inverse L^^ does not exist because of the zero 
eigenvalue of X3. This can be traced to the fact that the rotational U{1) symmetry ( unlike 
gauge symmetry ) can not be restored back to the original full SU (2) invariance. 
Thus we regularize L as follows 

2 4m^ 

^^^~N~ 'JP^ ^^^'^ ~ ^^^^^ 
Furthermore the value of in the matrix phase is expected to be very close to the classical 
value 1. This means in particular that 1 — ^ is a small number and thus it is a very good 
expansion parameter in this model. The vertex is given by 



iw? 1 11 
Vr, = ^(1 - + 4^(1 - ^)(1 - -)e.,3X3. (158) 



This is small both in expansion and in 1 — ^ expansion. We need to evaluate 
Tr^Tr^.i log [l + j^L-'V] = ^ ^-^_Tr;v-i(^- V)n..(^-'^)..^3--(^" V)...,. (159) 

r=l 
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Since 1) and V are (A^ — l)x(A^ — 1) ordinary matrices, 2) the trace jfziTr]\f_i becomes 
the ordinary integral in the hmit and 3) the operators Xa go over to the coordinates on the 
sphere we conclude that the only non-vanishing term in the — >oo limit is the order r = 1 
term in the above equation. We have 

TrsTr^v-ilogfl + ^L-V] = j^TrN_^{L-'V)i,i, + ... 

1 1 1 

= -4^(1 - -)eij3—TrN-iXi^ -—XjXs + ... 

(J) N x^ — 1 + e 

= 2" ,1 +-• (160) 

(p X3 — 1 + e 

This term is clearly going to zero in the limit. In particular the contribution of the zero 
eigenvalue of X3 is going to zero as ^^jz^- 

We can check that the quantum correction of the effective potential coming from the terms 
of order 0{jj) in equations (jlUjl - ipSj) are also going to zero in the limit. Hence the full quantum 
correction to the effective potential is given by pHjl . 



B The Star product on S| 

The coherent states on are constructed as follows. Let us introduce the 2— dimensional 
rank one projector Pi = ^ + Ua^. The requirement Pf = Pi implies the condition = 1. At 

the north pole we have uq = (0, 0, 1) and the projector becomes P? which projects onto the state 

2 

I 1 >= ^ Q ^ • III other words P? = |no, | >< rio, ||- A generic point n on is obtained 

by the rotation g such that n = gfiQ. The corresponding state is \n,^ >= g\nQ,^ > and the 

corresponding projector is precisely Pi which can also be rewritten as Pi = In, ^ >< n, ^|. 

2 2 ^ ^ 

The irreducible representation ^ of SU{2) can be obtained from the symmetric product of 
L copies of the fundamental representation |. The -I— representation of the element g E SU{2) 
is given by the matrix U'^^\g) defined by 

U^^\g) = g^s---®s9 , L times. (161) 

The (L + 1)— dimensional rank one projector Pl which defines the -I— coherent state |n, |^ > is 

given as the L— fold symmetric tensor product of the level \ projector Pi, viz 

^ 2 

Pl = \n, — >< n, — I = Pl (g)^...(g)s-Pi , L times. (162) 

The coherent state |n, ^ > can also be constructed as \n, ^ >= U'-^\g)\nQ, ^ > where |no, ^ > 
is the coherent state defined by the projector P£ = l^o, ■§ >< "^O; fl = P? (S>s...®s-Pi- 
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To any N x N matrix ( where = L + 1 ) we associate an ordinary function 0L(ri) on 
given by 

0i(n)=<n,^^|0|n,^>. (163) 
The product of two matrices (pip is mapped to the star product (pL *ipL{n) defined by 

(pL*'^L{n) =<n,^\(j)ip\n,^ > . (164) 

We can show that 

A:=0 1 ft i ft 

where 

Kab = Sab - naUb + ieabcflc- (166) 

Using these coherent states we can compute 

^ L, ,^ L L 

< n,-\La\n,- >= -Ha 

< n, ^\[La,(P]\n, I >= CaMn) (167) 



where Ca = —i^abcnbdc and 



^TrH= I ^<PL*Mn). (168) 



As another example we will compute -^TrLaAa which appears in the expansion (jH7|) of the 
one-plaquette action. We have immediately 



N -I f dVL, , , , , , , , , , 

--{na*{Aa)L+{Aa)L*na). (169) 



2N J 4:71 

It must be clear that {Aa)L is the function which corresponds to the N x N matrix Aa- In this 
article since we are mostly working with the N x N matrices Aa it is more easier to denote the 
corresponding functions {Aa)L ( in the very few places which appear ) by the same symbol Aa 
without fear of confusion . By using the star product p65j) we obtain the result 

§-JrLaAa = l^{<l>^ + l^daAa-^nada<!>^) 

= / ((1 + ^)'^- + l^daiAa - n.$oo)) . (170) 
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This is an exact formula where $00 is defined by $00 — f^a-^a- In the hmit $00 becomes exactly 
the normal component of and therefore Aa — Ua^oo is precisely the tangent gauge field on 
the sphere. Hence we can see directly that J dVLda{Aa — Ua^) = 0. Furthermore since $00 is a 
constant equal to in the limit we can conclude that we have the final result 

^TrL^A, = 0. (171) 

C The continuum limits of the one-plaquette actions Sp 

and S'p 

We need to check what happens to the action Sp in the large N limit, we have 

Sp = —Tr2{W + W+ -2). (172) 
A 

We will introduce the covariant matrices Aa and Sa defined respectively by 2Aa = Da — D'^ = 
Aa — A'^ and Sa = Da + D'a = 2{La + Aa) where Aa is a gauge field defined by the matrices 
2Aa — Aa-\- A^a- ^hc mcasurc becomes therefore dA'^dAa = 2^'^^ dAadAa- We start with the 
expansion 

7D = (^7 - ^l{^abc(^cFab) + ^l{eabc(^cFabf + 0{^)^ (173) 

and a similar expansion for 7^)/. We can now compute the first non- vanishing covariant terms 
in Tr2NW to be 

/ / 1 / / 1 / 

Tr2NW = Tr2Jv( 77 - -^77' (eabcf^c^ab) - J^^' ^^abc(^cFab) 

3 / 3 / / 

+ ^lli^abcCTcFabf + 1^11 {^abc(TcF'abf 

1 1 \ 

+ ^^l{eabcOcFab)l\eabcOcF'a,)+0{ — )y (174) 

Explicitly we have ( by reducing the 2A'"— dimensional trace Tr2N to the A^— dimensional trace 
Tr ) the following first contribution 

rr2jv77'+ h-c = ^Tr{l + ADaD'^)+ h.c 

32 

= ^N-—TrAl (175) 

Next we have 

Tr2N - ^11 {eabcCJcF'ab)^ + h.c = -^J'' (^abcDX, + eabcD'Xb + "^iFtaD'X^^ + h.c 

= -^,TT{^ea,cDcF'a, + ea,cDyab+{ADa,D[]-i[D,,D'S)F'^^ 

= -^JrTabF'a,. (17( 
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In above the matrices Tab are defined by Tab — Fab + F'ab ~ A^. Similarly we can obtain 

Tr2N{^- l{eabc(TcFa^ + ^.c = -^Tr^abFab- (177) 

Finally we need to evaluate the following three terms 

(3 / 3 / / 1 / / 

^J^J lieabcC^cFabf + 1^11 {^abcCrcF'aj^f + -J^j{eabcCrcFab)l {eabcC^cF'ab) 

= Tr2N (^-^l^ieabclJcFaby + l^li {^ahc(ycF'abf + -^{eabc(ycFoh)li {eabc(ycF'ab) 

+ -^[l,^ahc'7cFah]i{eahc(7cF'aS^. (178) 

We start by computing the last piece. To this end we use the identity 

[7, tabc(^cFab\ = ^(7a{^6, Fab] + ^CabcPc, Fab] = 2zaa{x^ , Fab} + O(^). (179) 

"0(;^)" stands for all other subleading terms which will yield corrections of the order of or 
higher to the action. The operators Xa are covariant coordinates on the fuzzy sphere defined 
by Xa — Ffalsfoi. It is clear that in the large N — >cxo limit x^ — >na which are the usual 
coordinates on the ordinary sphere. Thus the only difference between and the usual coor- 
dinates Xa = La/ \fc2 on the fuzzy sphere is that under U{N) gauge transformations we have 
Xa — ^gXaO^ as opposed to Xa which remain fixed. However since $ = the operator {xj^, Fab} 
tends in the continuum limit to 2nbFab which vanishes identically. Hence [7, eabcCcFab] = 
and thus we obtain 



(3 / 3 / / 1 / / 

^j;;^ll{eabcCrcFabf + l^ll {^abc(ycF'abf + J^{eabc(ycFab)li {eabcCycF'ab) 

+ (180) 
To evaluate the other terms we use the following remarkable identity 

7 7 = -^(D'aF>a + |) + -^eabcTc (^J'ab + l{Da, Ab} - l{Db, Aa} + i{D[, Ab} - i{Dl, 

= 4^(^-2Al-hAa,Sa]]+^eabcac(:Fa, + i{Sa,Ab}-i{Sb,Aa}). (181) 



or equivalently 



[7,7] = ^eabc<7c{Da, D',} + ^[Da, D'S (182) 
39 



and 



{7,7 } 



16 ^2 2 

2 - -^^a + J^eabcC^c-Fab- 



;i83) 



We see immediately that since we are already at order we can set in equation ()180|) the 
following 77'~1 and 7'7~1. Thus we obtain 

/ 3 3,2 /I 

5Sp + h.C = Tr2N ( J^ieabcC^cFabf + ^(eabcfTcFab)^ + ^(ea6cC^ci^a6)(eai,ct^ci^Ife) + 0{,J^, 



(184) 



The one-plaquette action Sp becomes ( by putting the contributions (jl75p . (jl76p . (jl77j) and 
(fTMj) together ) 



N 
J 



(185) 



We remark that -2TrTab{Fab + F'^,) +Tr{?>Fl + ^F'^, + 2FabF'^^) = Tr{Fab- F'J^ + %iTr{Fab + 
F'^^)[Aa,Ab] and thus 



N 
J 



- J^2TrAl + 4-Jr{Fab - F'J' + ^Jr{F^b + F:,)[A^, A] + 0(i 



;i86) 



By using the results Fab + F'^^ = ^[Sa,Sb] + 2i[Aa,Ab] + eabcSc and Fab - = i[Sa, Ab] + 
i[Aa, Sb] + 2eabcAc we have 



-^,Tr{Fab - F'a,f + ^Jr{Fab + F'a,)[Aa,Ab] 



Tr 2[Sa, Ab] [Sb, Aa] - 2[Sa, AbY + SA^ 



+ 16zeabcAc[Sa, Ab] - A[Sa, Sb][Aa, Ab] - lQ[Aa, Ab] 

(1 

We recall that Sa = 2La + 2Aa and that all commutators [A^, A^], [AajA^] and [AajA^] are of 
order and hence lead to terms of order in the action in the limit. With this approxima- 
tion the transformation laws Aa — >gAag^ and Aa — >gAag'^ + g[La,g^] become Aa — ^Aa and 
Aa — >Aa + g[La, g^] respectively. Thus we obtain 



-^Jr{Fab - F'J' + ^Tr{Fab + F'j[Aa,Ab] 



Tr 2[La, Ab] [Lb, Aa] - 2[L„ Ab]^ + 2AI 



16 



+ SieabcAc [La ,Ab]-A [La , Lb] [Aa ,Ab]) +0{ 



±). 

(188) 
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The one-plaquette action takes therefore the form 



N 



32 16 \ 1 

]^^^^a + i^t[La, A,] - t[L,, Aa] + eatcA^j + 0{— 



We find now the continuum hmit of the quadratic action 



S' 



N 



TTr2N{W^ + W+^ -2). 



(189) 



(190) 



We have 



where 



and 



3 I 3 ; / / / 

^llieabcC^cFabf + {eabc(rcF'^bf^ + 1 {^abc(ycF'^b)l{^abc(ycFab)- 



(191) 



(192) 



(193) 



Straightforward computation using equation (jl8ip gives 



Tr2^(7 7)' + h.c = ^Tr(^-4A3%^Tr[^,^J^- ATr(^{5.,A}-{^,,^}y 



4^ + (^^{Alf - N^Al^ + ^Tr S^f - ^Tr A,} - {5„ Aj) ' 

(194) 



Exphcitly we have 



Tr[{Sa.A,}-{Sb.A,} 



%Tr 



[Sa, A.f - [Sa, A,] [Sb, - 2[Sa, Sb] [A^, A] 



+ 4AlSl - AAbSbSaA, 
16 



- AAt 



Tr 



[Sa, A,Y - [Sa, A,] [Sb, - 2[Sa, Sb] A, 



+ ^JrAbSbSaA^ + ^Jr(^4{Air - N'Al 



(195) 
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In the last line above we have used the constraint = 4c2 — 4A^. By using the second 
constraint SaAa = —AaSa we can rewrite this equation as 



^Jr(^{S,,A,}-{S,,A,}^ = ^Tr 



[Sa, + [Sa, A] [Sb, Aa] + 2[Sa, S,] [A^, A,] 



+ ^Ai 



- i^Tr[^, + ^Tr {A{Alf - N^Al )(.196) 



Thus we obtain the final exact expression 

128„ ,o 512 



Tr2^(7 7)' + h-c = AN - ^-^TrAl+ ^-^Tr{Alf + ^Jr^l 

- [Sa. A,f + [Sa. A,][S,. Aa] + 2[5„, S,][Aa. A,] + 4^^ 



16^ 



(197) 



The next computation is to find 



(198) 



We use equation (jl81|l in the form ^'^ = I + eabcC^Jab = I + tabc(^c{^^ab + -^^Kab)- The 
definition of the operators /, Jab = —Jba and Kab = —K^a is of course obvious. Thus we can 
compute 



2 / 4 

- JpTr2N{l lY{,^abc(^cFab) + h.c = -^Treabc{Vc + V+)F„b. 

The operator Vc is defined in terms of / and Jab as follows 

Vc = ^ieabdJabJcd + ^abc{I J ab + JabI)- 



(199) 



(200) 



It is easy to check that the contribution of the first term 2itabdJabJcd is of order at least 
whereas the contribution of the second term eabdUab + JabI) is given by 



32 1 
—TrJ^^bFab + Oij^: 



The final result is 



2 32 1 

^Tr2^7V2 + h.c = -—TrJ^^,{Fab + F'j+Oi — 



(201) 



(202) 



Next we have to compute the following 

TFlTr2N7jh = TFlTr2N{^{l7)^{eabc(^cFabT + ^{llf{^abc(rcF'abf + {llf{eabc(ycFab){eabc(TcFc 



ab) 



+ ill {(^abcC^cF'ab) [7, eabcfycFab]^ 



(203) 
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Since [7, eabcO'cFab] is of order we obtain 



Ar4 



Tr2Nl ih 



3 ; 3 / / / / 

Tr2N[ -{.ilY{.(-abc(ycFabf + T^ill )i.^abcCr^Fl^f + {^^ f {eabcCrcKb)i.^abcCrcFab) 



(204) 



In above we can also make the approximations 77',7'7~1 since we are already at order 
Hence we obtain 



Ar4 • 



iV4 



Tr2Nl ih + h.c 



Tr2N K^^abcf^cFabf + "^ieabcC^ cF'ahY + ^^eabcC^cF'^b) i^^abcCTcFab) + 0{ 



—Tr2N ^Fl + ?,f2 + 2F,,F:, + 0( — 



(205) 



Finally we need to compute 

1 



1 



Tr2Nll + h-c 



Tr2N 7 7(eafecCrci^afe)7 li(^abc(^cFab) + 11 {.(^abc(^ cF^j,)ll (eabcCTcF^b) 



+ '^{liy{(^abc(^cFab){eabc(^cF'^b) + '^ll l{(^abc(^cFab)[l , ^abcC^ cF'ab]j + h-C 
= TFlTr2N ( {eabc(TcFabf + {eabc(TcF'^bf + 2{eabcCrcFab) (eabcCTcF'^b)^ + Oi^) 

(206) 



Ar4- 



^Tr2N{F!b + f2 + 2F„,F:,) + O(^). 



By putting equations (IT!I7|l . (IM) . (0nK|) and 

becomes 



together the quadratic one-plaquette action 



S' 



N 



128, 



512, 



16 



TrAi + j^Tr{AiY + — Tr ( - [^,, A,]^ + [^,, Ab][Sb. A,] + 4Af 



+ 2[^,, Sb][A^, Ab] ) + ^Jr{F^b - F'J' + ^Tr (F,, + F'J[A^, A] - ^Tr[^, A] 



iV4 



(207) 



As before if we drop all commutators [Aa, Ab], [Aa, Ab] and [Aa, Ab] ( since they are of order 
and hence lead to terms of order in the action ) then the limit of S'p reduces to 



S' 



256_ J ^ -jn bi'Z^ , J, ,j ^ t)4_ / - 

— TrA^L, + ^TrAl + j^Tr{Al + {A,, L J)^ + — Tr A, 



128, 



512, 



64 



(20^ 
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